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Resumen 


El diseño hidrológico de los embalses por construir o la revisión de los ya 
existentes requiere la estimación del llamado hidrograma de la creciente 
de diseño. La manera más simple y aproximada para estimar tal 


hidrograma, para un determinado periodo de retorno conjunto, es a 
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través del análisis de frecuencias bivariado del gasto máximo (Q) y el 
volumen escurrido (VW) anuales de las crecientes registradas. Las 
funciones Cópula (FC) son modelos probabilísticos basados en la 
dependencia entre Q y V, que establecen fácilmente su distribución 
bivariada con base en las funciones marginales previamente adoptadas o 
distribuciones de cualquier tipo, iguales o diferentes. La aplicación de las 
FC en las estimaciones hidrológicas mostró que un aspecto decisivo en su 
selección idónea está relacionado con la dependencia en el extremo 
derecho de los datos (A£'“) y la que tienen (Ay) ciertas FC. Por lo anterior, 
en este estudio se exponen las FC: t de Student, Gumbel-Hougaard, 
Clayton Asociada y Joe, que muestran valores de 1, que van en aumento. 
Se contrastan los valores de A, contra los 15“ obtenidos en 16 registros 
reales aleatorios de Q y V, para establecer la aplicabilidad de cada FC 
citada. Además, se procesa como aplicación numérica el registro de 26 
datos de Q y V anuales de las crecientes de entrada a la presa Adolfo 
López Mateos (Humaya) del estado de Sinaloa, México. Por último, se 
exponen las conclusiones, las cuales destacan las ventajas de las FC en 


los análisis de frecuencias bivariados de crecientes. 


Palabras clave: funciones Cópula, cociente tau de Kendall, coeficiente 
rho de Spearman, dependencia en el extremo superior, dependencia 


observada, periodos de retorno conjuntos, periodo de retorno secundario. 


Abstract 


The hydrological design of the reservoirs to be built or the revision of the 


existing ones, requires the estimation of the so-called Design Flood 
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Hydrograph. The simplest and most approximate way to estimate such a 
hydrograph, for a determined joint return period, is through the bivariate 
Frequency Analysis of the maximum flow (Q) and the annual runoff 
volume (V) of the registered floods. Copula Functions (FC) are 
probabilistic models based on the dependence between Q and V, which 
easily establish their bivariate distribution, based on previously adopted 
marginal functions or distributions of any type, equal or different. The 
application of the FC in the hydrological estimates showed that a decisive 
aspect in their ideal selection is related to the dependence on the extreme 
right of the data (25%) and that which have (2,) certain FC. Therefore, in 
this study the FC are exposed: Student's t, Gumbel-Hougaard, Clayton 
Associate and Joe, which show increasing values Ay. The values of 2, are 
contrasted against the 2“ obtained in 16 random real records of Q and 
V, to establish the applicability of each cited FC. In addition, the record of 
26 annual Q and V data of the inflow floods to the Adolfo López Mateos 
Dam (Humaya), in the state of Sinaloa, Mexico, is processed as a 
numerical application. Finally, the Conclusions are presented, which 
highlight the advantages of FC in the bivariate frequency analysis of 
floods. 


Keywords: Copula Functions, Kendall's tau quotient, Spearman's rho 
coefficient, dependence at the upper tail, observed dependence, joint 


return periods, secondary return period. 
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Introducción 


Generalidades 


La república mexicana se localiza en la zona de influencia de los huracanes 
que se originan en el mar Caribe y en el océano Pacífico. Además, lo 
afectan otros fenómenos meteorológicos de gran alcance espacial, como 
son los frentes fríos y las tormentas convectivas extensas, algunas de 
origen orográfico. Los fenómenos atmosféricos citados producen lluvias 
de gran magnitud, las cuales generan crecientes o avenidas máximas que 
inundan extensas regiones del país. Tales inundaciones causan pérdida 
de vidas humanas, y enormes daños económicos y ambientales. Este 
panorama que se repite cada año en distintas zonas de México destaca la 
importancia del estudio de las crecientes, con la finalidad de entender su 
dinámica e intentar cuantificar sus magnitudes, con el propósito de 
diseñar y proteger las obras hidráulicas de protección y aprovechamiento 
que se construyen (Aldama, 2000; Aldama, Ramírez, Aparicio, Mejía- 
Zermeño, 8: Ortega-Gil, 2006; Gómez, Aparicio, 8, Patiño, 2010). 


Las principales obras hidráulicas de protección son los diques y muros de 
contención que por lo general están relacionados con las rectificaciones y 
canalizaciones. Los puentes y otras obras del drenaje urbano también son 
estructuras de protección, cuyo diseño hidrológico se realiza con base en 


las crecientes de diseño (CD), que son gastos máximos del río asociados 
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con bajas probabilidades de ser excedidos. La estimación más confiable 
de las CD se realiza a través del análisis de frecuencias de crecientes 
(AFC), el cual intenta interpretar el registro de gastos máximos anuales 
disponible en términos de eventos futuros de una determinada 
probabilidad de ocurrencia (Bobée 8, Ashkar, 1991; Rao 8 Hamed, 2000; 
Meylan, Favre, 8 Musy, 2012). 


Para que los resultados del AFC sean aproximados y relativamente 
confiables se deben cumplir los siguientes cuatro requisitos: (1) el registro 
de gastos máximos anuales debe ser aleatorio; (2) la función de 
distribución de probabilidades (FDP) que lo represente debe ser idónea; 
(3) el ajuste de las FDP probadas se tiene que hacer con métodos 
estadísticos eficientes, y (4) la selección de resultados requiere ser 
objetiva (Hosking €: Wallis, 1997; Stedinger, 2017). 


Por otra parte, el diseño hidrológico de los embalses o presas de 
almacenamiento necesita la estimación del llamado hidrograma de la 
creciente de diseño, ya que en tales obras hidráulicas el volumen de la 
creciente es sumamente importante para establecer su peligrosidad y 
definir su manejo y almacenamiento temporal sin riesgo. Aldama (2000) 
demostró que los embalses no son sensibles al tiempo que trascurre en 
llegar al gasto máximo, pero sí a su magnitud y al volumen de la 
creciente. Lo anterior permite simplificar de forma notable el análisis 
multivariado de las crecientes, pues el enfoque bivariado del gasto 
máximo (Q) y del volumen (V) anuales permite definir con suficiente 
aproximación el hidrograma de la creciente de diseño, pues V tiene 


correlación positiva con la duración total (D) del hidrograma. 
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Los análisis bivariados de las crecientes comenzaron a finales del siglo 
pasado, con los trabajos de Goel, Seth y Chandra (1998), de Yue (1999), 
y de Yue, Ouarda, Bobée, Legendre y Bruneau (1999), que procesaron 
los datos conjuntos de Q y V, y de V y D. Al inicio de este siglo, Yue y 
Rasmussen (2002) establecen los conceptos hidrológicos básicos del 
análisis probabilístico bivariado. Los primeros análisis bivariados de 
crecientes emplearon la distribución Normal bivariada y el llamado modelo 
Logístico, que acepta funciones marginales de valores extremos, como las 
distribuciones Gumbel y GVE. Estos modelos tienen dos limitaciones 
importantes: (1) sus distribuciones marginales son iguales, y (2) los 
parámetros de ajuste de tales marginales intervienen en la modelación 
conjunta del registro de Q y V (Favre, El Adlouni, Perreault, Thiémonge, 
8, Bobée, 2004). 


En la actualidad, las limitaciones citadas son evitadas al aplicar las 
llamadas funciones Cópula (FC), modelos matemáticos que aceptan 
cualquier FDP como marginal, para formar la distribución bivariada o 
multivariada, y que se basan en la correlación que existe entre Q y V 
(Shiau, Wang, €: Tsai, 2006; Zhang € Singh, 2006). El uso de las FC se 
ha consolidado, como lo demuestran los textos especializados existentes 
(Salvadori, De Michele, Kottegoda, 8: Rosso, 2007; Zhang € Singh, 2019; 
Chen á Guo, 2019; Chowdhary 8 Singh, 2019). 


Ramírez-Orozco y Aldama (2000), Escalante-Sandoval y Reyes-Chávez 
(2002), así como Volpi y Fiori (2012), y Requena, Mediero y Garrote 
(2013) indican que el AFC bivariado conduce a una infinidad de 
combinaciones de gasto pico y volumen para una probabilidad de 


excedencia conjunta adoptada. Lo anterior implica que para un mismo 
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periodo de retorno conjunto existen muchas crecientes o hidrogramas que 
producirán distintos efectos en el embalse que se diseña o revisa, 
seleccionando lógicamente el que genera las condiciones más severas en 


su almacenamiento y vertedor de excedencias. 


Objetivos 


Los trabajos de Dupuis (2007), y de Poulin, Huard, Favre y Pugin (2007) 
demostraron que un aspecto decisivo en la selección de la FC idónea por 
utilizar está relacionado con la dependencia en el extremo derecho de los 
datos y la que muestran ciertas familias de FC. Por lo anterior, los 
objetivos básicos de este estudio fueron los cuatro siguientes: 1) exponer 
cuatro familias de FC que tienen dependencia significativa en su cola 
derecha (Au); 2) describir cómo se estima la Au del registro de Q y V 
disponible (Af*“); 3) contrastar los valores de Af“ obtenidos en 16 
registros reales de Q y V contra los que definen las cuatro FC expuestas, 
y 4) procesar como aplicación numérica el registro de 26 datos de Q y V 
anuales de las crecientes de entrada a la presa Adolfo López Mateos 
(Humaya) del estado de Sinaloa, México. 
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Funciones Cópula 


Concepto y definición 


Como ya se indicó, la ventaja esencial de las funciones Cópula (FC) 
consiste en permitir expresar una distribución conjunta de variables 
aleatorias correlacionadas, como una función de sus distribuciones 
marginales, previamente adoptadas. Entonces, una FC enlaza o relaciona 
las distribuciones marginales univariadas para formar una distribución 
multivariada. Otra ventaja básica de las FC al formar distribuciones 
multivariadas es el hecho de que ellas separan el efecto de la dependencia 
entre las variables aleatorias de los efectos de las distribuciones 


marginales en la modelación conjunta. 


Debido a lo anterior, la construcción de la distribución multivariada se 
reduce al estudio de la relación que guardan las variables correlacionadas 
si las distribuciones marginales univariadas son conocidas. El uso de las 
FC ofrece entera libertad para adoptar o seleccionar las distribuciones 
marginales univariadas que mejor representen a los datos (Shiau et al., 
2006; Kottegoda 8 Rosso, 2008; Meylan et a/., 2012; Zhang € Singh, 
2019). 


Como en este estudio aplicarán las FC en el análisis de frecuencias 
bivariado de las crecientes anuales, la definición siguiente se refiere a dos 
variables aleatorias X y Y correlacionadas, cuya función de distribución de 


probabilidades acumuladas conjuntas es Fxyv(x,y), con distribuciones 
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marginales univariadas Fx(x) y FMy); entonces la FC existe, es C[:] y es 


tal que: 
FxyGco y) = CIF), Fy(y)] (1) 


La ecuación anterior define el concepto básico para el desarrollo de las FC 
y se conoce como Teorema de Sklar, expuesto en 1959 (Shiau et al., 
2006; Kottegoda € Rosso, 2008; Meylan et a/., 2012; Zhang € Singh, 
2019). 


Familias de Cópulas 


Las funciones Cópula (FC) que han sido desarrolladas se han clasificado 
en cuatro clases: de Arquímedes, de valores extremos, elípticas y 
misceláneas. También se clasifican en Cópulas de un parámetro o de 
varios, dependiendo de la amplitud con la cual la estructura de la 
dependencia entre las variables Q y V es definida (Meylan et a/., 2012; 
Chowdhary 8 Singh, 2019). 


Salvadori et al. (2007) exponen un resumen amplio y útil de FC que han 
sido aplicadas en el campo de la hidrología. Al respecto, en la relación 


siguiente no se han incluido las familias de Cópulas de Ali-Mikhail-Haq, 
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llamada Cópula AMH, debido a que limita la dependencia a un rango del 
cociente tau de Kendall de -0.1817 a 0.3333, y la familia de Farlie- 
Gumbel-Morgenstern, designada FGM, que lo limita aún más, al intervalo 
de +0.2222. 


Las Cópulas de Arquímedes han tenido aplicación amplia debido a su 
construcción simple, un único parámetro, rango amplio y aceptación de 
ambos tipos de dependencia (positiva y negativa). Designando a FxQO)=u, 
FAy)=v y 0 el parámetro que mide la dependencia o asociación entre u y 
v, se tienen las siguientes dos familias de Cópulas de Arquímedes, las 
cuales aceptan dependencia negativa y positiva (Shiau et al., 2006; 
Genest 8 Favre, 2007; Salvadori et al., 2007; Michiels £ De Schepper, 
2008; Zhang 8 Singh, 2019; Chen € Guo, 2019; Chowdhary € Singh, 
2019). 


1. Clayton. Esta función Cópula es llamada de Cook-Johnson por Zhang 
y Singh (2006). Sraj, Bezak y Brilly (2015) la citan con los dos nombres. 


Su ecuación y espacio de variación de 0 son: 
_ (9 “ib -1/0 
C(u,v) = (u +v*?- 1) [—1, 00) M £03 (2) 


Para la dependencia positiva 9 > O y para la negativa -1<86<0, con O = 
O para la independencia entre u y v. La relación de 6 con el cociente tau 


de Kendall es la siguiente: 
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Tan = 515 (3) 


2. Frank. Su ecuación y espacio de variación de O son: 


(tua) (0 


Cu) = -5 In [1 + EE 25) (—oo, 00) N (0) (4) 


Para la dependencia negativa 0< 0 < 1 y para la positiva 9 > 1,conO = 1 
para la independencia entre u y v. La relación de 6 con el 7, es la 


siguiente: 
tn = 1+5[D,(8) - 1] (5) 
siendo D:(0) la función Debye de orden 1, cuya expresión es: 


D,(0) =3f, ds (6) 


es-1 


La ecuación anterior se estimó con integración numérica, ratificando sus 


resultados con los valores tabulados por Stegun (1972). 


Enseguida se cita una familia de Cópulas de valores extremos, la cual 


acepta sólo dependencia positiva. 
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3. Gumbel-Hougaard. Su ecuación y espacio de variación de 6 son: 
C(u,v) = exp [E Inu)? + (In Dye, [1, 00) (7) 


Con 6 = 1 se tiene independencia entre u y v. La relación de 6 con el 


cociente tau de Kendall es la siguiente: 


0-1 
e AT (8) 
Por último, se citan dos familias de Cópulas de clase miscelánea. 


4. Plackett acepta dependencia negativa (0 < 1) y positiva (9 > 1), con 


8 = 1 para la independencia entre u y v: 


14(0-1)(utv)  y[1+(0-1)(u+v)]?-4uvo(0-1) 


Cu, v) = 2(0-1) 2(0-1) 


con0<86<w (9) 


La relación de 6 con la rho de Spearman es la siguiente: 


_ (6+1)  261n(0) 
Pn= y e-1? 


(10) 
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5. Raftery acepta sólo dependencia positiva, y O varía entre O y 1, con 


8 = 0 para la independencia entre u y v: 


C(u,v) = min(u, v) + (wr — [max(u, v)]4+90/0-9) (11) 


Las relaciones de 0 con el cociente tau de Kendall y con el coeficiente rho 


de Spearman son las siguientes: 


n=3 (12) 
0(4-30) 
n —_ (2-0)? (13) 


Las familias de las Cópulas elípticas son distribuciones y sus marginales 
que aceptan la propiedad de ser elípticas, como las distribuciones Normal, 
t de Student y Cauchy. Un vector aleatorio X e R* se dice que tiene una 
distribución elíptica si admite la siguiente representación estocástica: X = 
u +RAU; en la cual, He RY, R es una variable aleatoria independiente de 
U, que es un vector aleatorio uniforme en la esfera unitaria en RY y A es 
una matriz fija de dxd, tal que AA” es no singular (Salvadori et a/., 2007; 
Chen 8 Guo, 2019). 


Cuando las distribuciones marginales son de tipos diferentes y 
posiblemente no elípticas, las Cópulas elípticas se denominan Cópulas 


Metaelípticas, cuyas familias principales son la Gaussiana y la de la t de 
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Student (Chowdhary 8, Singh, 2019). La Cópula t de Student será descrita 


con detalle posteriormente. 


Familias de Cópulas asociadas 


Las cuatro clases de familias de Cópulas que existen y para el caso 
bivariado tienen tres familias asociadas, se definen por las siguientes 
transformaciones (Salvadori €: De Michele, 2004; Michiels 8. De Schepper, 
2008; Chowdhary 8 Singh, 2019): 


C'(u,v) =u-—C(u,1-—v) (14) 
C"(u,v) =v-C(1—u,v) (15) 
É(u,v)=u+v-1+C(1-u 1-0) (16) 


Las dos primeras transformaciones cambian la naturaleza de la 
dependencia de negativa a positiva y viceversa. La tercera transformación 
se denomina Cópula de supervivencia (Survival Copula), por estar 
relacionada con la distribución conjunta de las probabilidades de 
excedencia, del par de variables aleatorias u y v cuya Cópula es C. 
Además, las tres familias de Cópulas asociadas tienen el mismo valor del 


cociente tau de Kendall; las dos primeras con su signo cambiado. 
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Dupuis (2007) destacó la importancia de la Cópula de Clayton Asociada 
por tener dependencia significativa en su cola derecha definida por la 
Ecuación (17), y cuya fórmula de acuerdo con las ecuaciones (2) y (16) 


es la Ecuación (18): 


Au = golpe (17) 


-1/0 


C(uv) =u+v-1+[(1-u)9+(1-v)2-1] (18) 


Integración numérica 


Para cuantificar varias integrales, por ejemplo, la Ecuación (6), se realizó 
una integración numérica con base en el método de la cuadratura Gauss- 
Legendre, cuya ecuación operativa es (Nieves € Domínguez, 1998; 
Campos-Aranda, 2003): 


IA (19) 


en la cual w; son los coeficientes del método cuyas abscisas son h; y np el 
número de parejas en donde se evalúa la función f(x), con el argumento 


indicado en f [-]. En Davis y Polonsky (1972) se obtuvieron las 12 parejas 
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utilizadas de w; y h; con 15 dígitos, que son aceptables en el leguaje Basic 


como variables de doble precisión. 


La ecuación operativa de la integración bivariada es la integral doble 


siguiente (Nieves 8 Domínguez, 1998): 


(b-aY(d-c) (b-a)h¡+b+a (d—c)h¡+c+d 


IC E 


En Davis y Polonsky (1972) se obtuvieron las 20 parejas utilizadas en esta 


integral doble de w; y hi con 15 dígitos. 


Medidas de asociación 


Concordancia 


Como la FC caracteriza la dependencia entre las variables aleatorias u y 
v, es necesario el estudio de las medidas de asociación para disponer de 
un método que permita estimar su parámetro 6. En términos generales, 
una variable aleatoria es concordante con otra cuando sus grandes 
valores están asociados con los grandes valores de la otra y los valores 
pequeños de una con los valores reducidos de la otra (Salvadori et al., 
2007; Chowdhary 8 Singh, 2019). 


Unas variables con correlación lineal directa serán concordantes, pues al 


aumentar una la otra también lo hace. Variables con correlación lineal 
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inversa serán discordantes, ya que a grandes valores de una le 
corresponderán pequeños valores de la otra. Entonces, se dice que (xi, 
y1) y QQ, y2) son concordantes si x1 < X2 Y Y1< Y2, O bien X1 > X2 y Y1 > 
y2. Son discordantes si Xi < X2 y Y1 > Y2, O bien x1 > X2 Y y1i< y2. Lo 
anterior implica que son concordantes (c) las parejas (x1-X2) (y1-Y2)>0 y 
discordantes (d) cuando (x1-X2) (yi1-y2)<0 (Salvadori et al., 2007; 
Chowdhary 8 Singh, 2019). 


Una medida numérica de asociación es un estadístico que indica el grado 
de dependencia o asociación de las variables. Para fines de comparación, 
tales medidas varían de O a +1 0 a -1, indicado la asociación perfecta 
positiva en +1 o negativa en -1. El cociente tau de Kendall y el coeficiente 
rho de Spearman son dos medidas no paramétricas que proporcionan 
información sobre una forma especial de asociación o dependencia, 


conocida como concordancia (Salvadori et al., 2007; Chen € Guo, 2019). 


Cociente tau de Kendall 


Mide la probabilidad de tener parejas concordantes, por lo cual es el 
cociente de c-d entre c+d. Su expresión para estimarlo con datos 
bivariados es (Zhang €: Singh, 2006; Zhang €: Singh, 2019): 


n= Ei signol( — xy) (0 — 9,)] (21) 


136 


Tecnología y ciencias del agua, ISSN 2007-2422, 


Open Access bajo la licencia CC  BY-NC-SA 4.0 | 14(5), 120-188. DOI: 10.24850/j-tyca-14-05-03 


(https: //creativecommons.org/licenses/by-nc-sa/4.0/) 


a 0) Check for updates 
OPEN ACCESS 


Tecnología y 


CienciaszAgua 


En la ecuación anterior, n es el número de observaciones y el signo[-] es 


+1 si tales parejas son concordantes y -1 si son discordantes. 


Genest y Favre (2007) exponen una prueba para el cociente tau, 
aceptando que la hipótesis nula Ho indica que X y Y son independientes y 
entonces el estadístico tiene distribución aproximadamente Normal con 
media 0 y varianza 2(2n+5)/[9n(n-1)]. Entonces, Ho será rechazada con 


un nivel de confianza a = 5 % si: 


9n(n-1) 
2(2n+5) 


Itn| > Zajz = 1.96 (22) 


Coeficiente rho de Spearman 


Mide la correlación entre parejas de rangos (R;, Si) de las variables 
aleatorias X; y Y¡. Por ello equivale al coeficiente de correlación de Pearson. 
Su expresión para estimarlo en un registro bivariado de tamaño n es la 
siguiente (Chowdhary 8 Singh, 2019; Zhang 8 Singh, 2019): 


Z Ri (23) 


nan ta e 


Salvadori et al. (2007) presentan otra fórmula para su cálculo a partir de 
un registro bivariado con n datos, la cual coincide con la utilizada por la 
WMO (1971): 
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Genest y Favre (2007) también presentan una prueba similar para el 
coeficiente rho, cuya distribución se aproxima a la Normal con media O y 


varianza 1/(n-1); por lo cual se rechaza Ho con un nivel a = 5 %, si: 


vn — 1-|p,| > 1.96 (25) 


Estimación del parámetro de dependencia 


El método más simple para estimar el parámetro 6 de las FC se asemeja 
al método de momentos y se basa en la inversión de la ecuación que 
relaciona a 6 con el cociente tau de Kendall o con el coeficiente rho de 
Spearman. Los otros dos procedimientos disponibles se denominan: 
método de máxima pseudo-verosimilitud y método exacto de máxima 
verosimilitud (Meylan et a/., 2012; Chowdhary € Singh, 2019; Zhang € 
Singh, 2019; Chen 8 Guo, 2019). Para obtener 6, en las ecuaciones (5), 
(10) y (13) se procede por tanteos; en cambio, en las ecuaciones (3), (8) 


y (12) se despeja su valor. 
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Estimación de probabilidades empíricas conjuntas 


Las probabilidades empíricas bivariadas se estimaron con base en la 
fórmula de Gringorten, aplicada por Yue (2000b), Zhang y Singh (2019), 
y Chen y Guo (2019). Tal fórmula es: 


10.44 (26) 


— n+0.12 


En la cual, ¡es el número de cada dato cuando están ordenados de manera 
progresiva, y n es el número total de ellos o amplitud del registro 
procesado. La expresión anterior se aplicó en el plano bidimensional, con 
los datos ordenados en forma progresiva; los gastos máximos (Q) en los 
renglones y los volúmenes (V) en las columnas. El plano formado es un 
cuadrado de n por n casillas, con n casillas en su diagonal principal, 
cuando el número de orden del renglón es igual al de la columna. 
Después, cada pareja de datos anual (Q y V) se localiza en el plano 
bidimensional, y la casilla definida por la intersección del renglón y 
columna se identifica con el número /, que corresponde al año histórico 
dibujado. 


Cuando las n parejas de datos están dibujadas, se busca el año 1 y se 
define un área rectangular o cuadrada de valores menores de Q y de V, 
cuyo conteo de casillas numeradas dentro es NM o combinaciones de Q 


y V menores. Calculados los n valores de NM;, se aplica la fórmula de 
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posición gráfica de Gringorten para estimar la probabilidad empírica 


conjunta o bivariada: 


F(xy)=P(0<qV<v) = 22 (27) 


n+0.12 


Selección de la Función Cópula 


Un enfoque simple de selección de la función Cópula se basa en los 
estadísticos del error de ajuste, al comparar las probabilidades empíricas 
observadas (w?) con las teóricas calculadas (wf) con la FC que se prueba. 
Los indicadores aplicados son el error medio estándar (EME), el error 
medio absoluto (EMA) y el error absoluto máximo (EAM); sus expresiones 
son (Chowdhary € Singh, 2019): 


EME = B n (w? —w£)? (28) 


EMA = ¿Eiza 1? — wp)! (29) 
EAM = max] (wp = w¿)| (30) 
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Dependencia en el extremo superior 


Generalidades 


El criterio más importante que se aplica para seleccionar una FC es el 
basado en la magnitud de la dependencia en el extremo superior de la 
distribución conjunta, lo cual tiene impacto en la veracidad de las 
predicciones extremas. La dependencia en la cola superior derecha (A) 
es la probabilidad condicional de que Y sea mayor que un cierto percentil 
(s) de Fy(y), dado que X es mayor que tal percentil en Fx(x), conforme s 
se aproxima a la unidad. La dependencia en la cola inferior izquierda (A, ) 
compara que Y sea menor que X, cuando s se aproxima a O (Chowdhary 
8: Singh, 2019; Salvadori et al., 2007). 


En relación con las FC expuestas, las de Frank y Plackett tienen 
dependencias insignificantes en sus zonas extremas, por ello A, = O y Au 
= 0. Por el contrario, la Cópula Gumbel-Hougaard tiene dependencia 


significativa en la cola superior igual a: 
ly =2- 219 (31) 
Las Cópulas de Clayton y Raftery la tienen en su cola inferior e igual a: 


(32) 
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(33) 


Como ya se indicó, la Cópula de Clayton Asociada tiene Au = 2/8 
(ecuación (17)). Familias de Cópulas no expuestas, con A, > 0 son las de 
Galambos, Genest-Ghoudi y Húsler-Reiss (Salvadori et al., 2007; 
Chowdhary 8 Singh, 2019). 


Dupuis (2007) probó seis familias de Cópulas y encontró que su capacidad 
para estimar eventos extremos varía de mala a buena, con el siguiente 
orden: Clayton, Frank, Normal, t-Student, Gumbel-Hougaard y de Clayton 
Asociada (Survival Clayton). A conclusiones similares llegan Poulin et al. 
(2007), al comparar las mismas seis familias de Cópulas y la denominada 


A12, que tiene dependencia significativa en su cola derecha. 


Estimación de la dependencia observada 


Para abordar la estimación de la dependencia en la cola superior (Au) que 
muestran los datos disponibles, se debe definir primero la llamada Cópula 
Empírica. Como la FC caracteriza la dependencia entre las variables 
aleatorias X y Y, entonces el par de rangos Ri; y S; procedentes de las tales 
variables son el estadístico, que retiene la mayor cantidad de información 
y su escalamiento con el factor 1/(n+1) genera una serie de puntos en el 
cuadrado unitario [O, 1]?, formando el dominio de la Cópula Empírica 
(Chowdhary €: Singh, 2019), definida así: 
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1 Ri Si 
Ca kt v) = a del E < el < v) (34) 


En la ecuación anterior, 1(+) indica una función de las variables aleatorias 
U y V, que son una transformación en aumento continuo de X y Y en 
relación con las integrales de probabilidad empírica FAX) y Fa(Y), cuyas 


ecuaciones son: 


__ Rango (X¡) 
TAR 


Rango (Y;) 


U; 
n+1 


FE, (xi) V¡ = = F, (Y) (35) 


Poulin et al. (2007) utilizan el estimador propuesto por Frahm, Junker y 
Schmidt (2005), el cual se basa en una muestra aleatoria obtenida de la 


Cópula empírica; su expresión es: 


1 E AL l 
dl =2- Z2exp e A In NO Gimnl (36) 


Este estimador acepta que la FC se puede aproximar por una de la clase 
de valores extremos; tiene la ventaja de no requerir un valor umbral para 
su estimación, como los cuatro expuestos por AghaKouchak, Sellars y 
Sorooshian (2013). 
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Familias de Cópulas por contrastar 


Justificación 


Desde los inicios de las aplicaciones de las funciones Cópula (FC) en los 
análisis hidrológicos multivariados, Favre et al. (2004) formularon en sus 
conclusiones que la etapa crucial de tales análisis radica en la selección 
adecuada de la FC a los datos disponibles. Posteriormente, Requena et al. 
(2013) destacaron que una FC que muestra un buen ajuste al total de 
datos disponibles no garantiza una representación adecuada de los 
valores extremos de tales datos y, por ello, hay que estimar la 
dependencia observada A£F“ (ecuación (36)), para buscar una FC que la 


reproduzca. 


Por lo anterior, cuando se requieren estimaciones de periodos de retorno 
conjuntos elevados, superiores a tres veces la amplitud del registro 
conjunto disponible de gastos y volúmenes máximos anuales, se debe 
buscar y seleccionar la FC que mejor reproduce la dependencia en la cola 
derecha de los datos (A5"“). En este estudio se contrastan cuatro familias 
de Cópulas, cuya dependencia en su extremo derecho varía de mayor a 
menor grado, las cuales son las siguientes: Joe, Clayton Asociada, 
Gumbel-Hougaard y t de Student. Las tres primeras pertenecen a la clase 


de las Cópulas de Arquímedes y la última es de la clase elíptica. 
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Cópula de Joe 


La ecuación básica de esta familia de FC, que sólo acepta dependencia 
positiva, es la siguiente (Joe, 1993; Salvadori et al., 2007; Chowdhary 8 
Singh, 2019): 


1/0 


C(u,v) =1-|(1-w+(1-v- (1-w0%41- vw] (37) 


cuyo parámetro de dependencia es O > 1, con O = 1, para el caso de 
independencia entre u y v. Como no tiene una expresión que relacione su 
parámetro 6 con el cociente tau de Kendall, éste se estima en función de 
su ecuación del generador ¿(s) y de su derivada $'(s) de la FC, cuya 
expresión es (Salvadori et al., 2007; Chowdhary 8, Singh, 2019): 


19(5) 
TA=1+4f, amo ds (38) 


en la cual s es la variable aleatoria unitaria 0<s < 1. Las ecuaciones de 
su generador y de la derivada de éste son (Salvadori et a/., 2007; Zhang 
8: Singh, 2019): 


p(s) = —In[1 — (1 — sy?] (39) 
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9(1=5)9-1 


PS) == a. (40) 


Teniendo como dato el valor de cociente tau de Kendall, la Ecuación (38) 
se integró numéricamente con base en la Ecuación (19), para obtener por 
tanteos el valor de O que la cumple. Obtenido el valor de 6, la dependencia 
significativa en el extremo derecho (A) de la Cópula de Joe, se obtiene 


por medio de la Ecuación (31) de la Cópula de Gumbel-Hougaard. 


Cópula de Clayton Asociada 


La Cópula definida con la Ecuación (18) será aplicada y contrastada 
debido a que Dupuis (2007) la expuso y mostró que tiene una fuerte 


dependencia significativa en su extremo derecho (2). 


Cópula t de Student 


Demarta y McNeil (2005) indican que las Cópulas de Gumbel-Hougaard y 
Clayton Asociada son de aplicación más simple que las Cópulas extremas 
derivadas de la t-Student, pero esta última define un valor inferior de la 
dependencia, como lo señaló previamente Dupuis (2007) y lo ratificó 
Poulin et al. (2007). La ecuación de la Cópula t de Student es la siguiente 
(Salvadori et al., 2007; Chowdhary 8 Singh, 2019): 
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tuu) ptit(v) 1 x2-20:x"y+y? -(v+2)/2 
C(u, v) = 7 NA El A) dx: dy (41) 


En la expresión anterior v > 2, O varía entre -1 y 1, con O = O para la 
independencia entre u y v; t;*(:) es la inversa de la distribución t de 
Student univariada, con v grados de libertad, para la probabilidad indicada 


en el paréntesis. 


Entonces, con la idea de acotar la dependencia en la cola derecha en un 
valor inferior se exponen las ecuaciones de la Cópula t-Student (Demarta 
8: McNeil, 2005; Poulin et a/., 2007; Salvadori et al., 2007; Genest, Favre, 
Béliveau, 8i Jacques, 2007): 


Tn =“arcsen 0 (42) 


T 


VVFIVIO 
ho 2 (E) (43 


El signo menos dentro del paréntesis indica que tal valor es la probabilidad 
de excedencia (survival probability) que debe definir la variable aleatoria 
con distribución t de Student y v + 1 grados de libertad. Para evitar el uso 
de algoritmos numéricos que estimen el valor de t,.,, en función de la 
probabilidad de excedencia, Demarta y McNeil (2005) exponen la 


siguiente fórmula alternativa: 
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En la ecuación anterior, Á no lleva subíndice debido a que la Cópula t de 
Student tiene igual dependencia significativa en sus dos colas o extremos. 
La Ecuación (44) se integró numéricamente con base en la Ecuación (19) 


y con tales resultados se elaboró la Tabla 1, que es una ampliación de la 
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7/2 
eN _ des 0)/2 
1/2 


o 


cos's:ds 


cosYs-ds 


OPEN Que: 


expuesta por Demarta y McNeil (2005). 


0) Check for updates 


Student para diversos valores del y y 0. 


(44) 


Tabla 1. Valor de la dependencia en los extremos de la Cópula t de 


Valor del parámetro de asociación (0) 


-0.40 


-0.50 


-0.75 


0.00 


0.40 


0.50 


0.55 


0.60 


0.65 


0.70 


0.75 


0.80 


0.85 


0.90 


0.95 


0.98 


0.0773 


0.0577 


0.0195 


0.1817 


0.3393 


0.3910 


0.4195 


0.4502 


0.4833 


0.5195 


0.5594 


0.6042 


0.6557 


0.7177 


0.7995 


0.8729 


0.0378 


0.0189 


0.0257 


0.0117 


0.0061 


0.0020 


0.1161 


0.0756 


0.2606 


0.2031 


0.3125 


0.2532 


0.3419 


0.2822 


0.3739 


0.3144 


0.4091 


0.3503 


0.4481 


0.3907 


0.4918 


0.4366 


0.5415 


0.4896 


0.5995 


0.5523 


0.6702 


0.6298 


0.7648 


0.7349 


0.8505 


0.8311 


0.0049 


0.0025 


0.0002 


0.0331 


0.1269 


0.1705 


0.1970 


0.2275 


0.2625 


0.3031 


0.3506 


0.4071 


0.4758 


0.5630 


0.6845 


0.7980 


0.0013 


0.0006 


0.0000 


0.0150 


0.0811 


0.1173 


0.1405 


0.1679 


0.2004 


0.2393 


0.2861 


0.3434 


0.4151 


0.5086 


0.6424 


0.7699 


10 


0.0004 


0.0001 


0.0000 


0.0069 


0.0527 


0.0819 


0.1015 


0.1255 


0.1549 


0.1911 


0.2360 


0.2925 


0.3652 


0.4627 


0.6059 


0.7452 


y = grados de libertad. 


En la Tabla 1 se destaca que incluso para valores negativos y correlación 


O, la Cópula t de Student tiene dependencia asintótica en sus extremos y 


que tal dependencia decrece con el aumento de los grados de libertad. 
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La aplicación o integración de la Cópula t de Student (Ecuación (41)) se 
realizó con base en la Ecuación (20), con 20 parejas de coeficientes w; y 
abscisas h; de 15 dígitos tomados de Davis y Polonsky (1972). Para la 
estimación de los límites superiores de las integrales se aplicó el algoritmo 
expuesto por Zelen y Severo (1972) para estimar la variable aleatoria t 
de Student con v grados de libertad, en función de la variable Normal xp, 


en la cual p es la probabilidad de excedencia: 


box + OD y BAD y O y 20) (45) 
siendo: 

90) =(7+x) (46) 
g2(x) = (5x5 + 16x?* + 3x) (47) 
alt= Gx + 19x* + 17x? — 15x) (48) 
94(%) = 5 (091? + 776x7 + 1482x5 — 19201? — 945x) (49) 


Cuando p < 0.50 hacer t = -t. El algoritmo para estimación de la variable 


Normal xp también procede de Zelen y Severo (1972) y es el siguiente: 
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Xp =0 


en la cual: 

= pl 
¿= ME 
siendo: 
Co = 2.515517 
c1 = 0.802853 
c2 =0.010328 
dí = 1.432788 
dz = 0.189269 
d3 =0.001308 


Sip > 0.50 hacer p = 1 - p y Xp = -Xp. Los algoritmos definidos por las 


ecuaciones (45) a (51) fueron probados con los valores numéricos que 


Co+ci:9+c2:0? 
1+d1:8+d,:92+d3:83 


OPEN Que: 


exponen Zhang y Singh (2019). 


Para establecer los límites inferiores de las integrales de la Ecuación (41) 
se utilizó la tabulación de la distribución t de Student, expuesta por Ostle 
y Mensing (1975) como Apéndice 5. En la Tabla 1 se definen para una 


probabilidad de no excedencia de 0.9995 los siguientes valores de t con 


Open Access bajo 
(https: //creativecommons.org/licenses/by-nc-sa/4.0/) 


la 


licencia 


cc 


BY-NC-SA 


4.0 


0) Check for updates 


(50) 


(51) 
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vgrados de libertad de 2 a 10: 31.598, 12.941, 8.610, 6.859, 5.959, 
5.405, 5.041, 4.781 y 4.587. Tales valores con signo menos se tomaron 
como límites inferiores según el grado de libertad asignado en la 
aplicación de la Ecuación (41), el cual, lógicamente, varió de 2 a 10; 
seleccionado la Cópula t de Student de mejor ajuste (ecuaciones (28) a 


(30)). 


Ratificación de la función Cópula seleccionada 


Esta es la etapa más importante del proceso de aplicación práctica de las 
FC, pues en ella se verifica que tal modelo reproduzca fielmente las 
probabilidades conjuntas observadas (Ecuación (27)). Yue (2000a) indica 
que la forma simple y práctica de representar las probabilidades conjuntas 
empíricas y teóricas consiste en llevar al eje de las abscisas la primera y 
al eje de las ordenadas la segunda; lógicamente, en tal gráfica, cada 
pareja de datos define un punto que coincide o se aleja de la recta a 450, 
La inspección de la gráfica descrita y el valor del coeficiente de 
correlación, en estos casos, superior a 0.98, ratifican la validez del modelo 


probabilístico conjunto. 


Yue (2000b), y Yue y Rasmussen (2002) aplican el test de Kolmogorov- 
Smirnov, con un nivel de significancia (a) del 5 %, para aceptar o rechazar 
la diferencia máxima absoluta (dma) entre las probabilidades conjuntas. 
Para evaluar la estadística (D,) del test se utilizó la expresión que exponen 
Meylan et al. (2012), paraa = 5%: 
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Dn == (52) 


n es el número de datos. Si la dma es menor que D, se ratifica la FC 


adoptada. 


Conceptos del análisis de frecuencias bivariado 


Periodos de retorno bivariados 


El primer periodo de retorno bivariado del evento (X, Y) se define bajo la 
condición OR, lo cual indica que los límites x o y, o ambos pueden ser 
excedidos, y entonces la ecuación clásica del periodo de retorno o inverso 
de la probabilidad de excedencia será (Shiau et al., 2006; Genest 8 
Chebana, 2017): 


a ———_——— A (53) 


T PQoxorY>y) I-Fxy(0y)  1-C[Fx()Fy0)] 


en la cual C[Fx (Go), Fy(y)] es la FC seleccionada. 


El segundo periodo de retorno bivariado del evento (X, Y) está asociado 
con el caso en que ambos límites son excedidos (X > x,Y > y) o condición 
AND; su ecuación es (Shiau et a/., 2006; Genest 8: Chebana, 2017): 
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(54) 


xr oa ran EEN RO AO CERDO] 


Aldama (2000) obtiene la expresión Fx y(x, y) de la probabilidad bivariada 
de excedencia mediante un razonamiento de probabilidades lógico y 
simple aplicado en el plano cartesiano. En cambio, Yue y Rasmussen 
(2002) recurren al plano cartesiano para definir el evento bivariado (X, Y) 


conceptualmente, que puede ocurrir en alguno de los cuatro cuadrantes. 


La relación entre los periodos de retorno bivariados y los univariados es 
la siguiente (Yue € Rasmussen, 2002; Shiau et al., 2006; Vogel 8 
Castellarin, 2017): 


Txy < min[Tx, Ty] < maxI[Tx, Ty] < T'xy (55) 
siendo: 
o 6 > 6 NA0) 


Salvadori y De Michele (2004) introducen con detalle el concepto del 


periodo de retorno secundario (7), designado así para enfatizar que el 
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periodo de retorno conjunto Txy es el primario, de donde procede al utilizar 


las isolíneas que define la FC aplicada, cuya expresión es: 
Ls = [((u,v) € F:C(u,v) = 5] (58) 


siendo s la variable aleatoria unitaria 0 <s < 1 y C la FC probada. Después 
se define una región Bc(s) en el espacio unitario (1?) sobre la isolínea, 


debajo de ella y a la izquierda, que será: 
B¿(s) = [(u,v) € 12: C(u,v) < s) (59) 


En FC de la clase de Arquímedes, la distribución univariada de Kendall, 
designada Kc(s), proporciona una medida de los eventos dentro de la 
Bc(s); su ecuación es (Salvadori 8 De Michele, 2004; Salvadori 8 De 
Michele, 2007; Salvadori et al., 2007; Gráler et al., 2013): 


K¿(s) =s-= 2 (60) 


en la cual p(s) es el generador de la FC y 9'(s) su derivada. Por último, el 


periodo de retorno secundario (7) de los eventos fuera de Bc(s) es: 


¿== (61) 


— 1EKC(S) 
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cuyo denominador es la probabilidad de excedencia (survival function), la 
cual corresponde a eventos posiblemente destructivos o peligrosos. La 
relación que guarda Y con los periodos de retorno conjuntos es la siguiente 
(Gráler et a/., 2013): 


Txy €S STxy (62) 


La aplicación de la Ecuación (61) requiere la estimación previa de la 
Ecuación (60) para las siguientes seis probabilidades de no excedencia 
(s): 0.98, 0.99, 0.998, 0.999, 0.9998 y 0.9999, que corresponden a los 
periodos de retorno conjuntos de 50, 100, 500, 1 000, 5 000 y 10 000 
años. En las Cópulas de Gumbel-Hougaard y de Joe, las ecuaciones de su 


generador y su derivada son las siguientes (Zhang y Singh, 2019): 


p(s) = Ins) (63) 
$5 =-¿Einsyen (64) 
p(s) = —In[1 — (1 — sy?] (65) 
9 = (66) 
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Eventos críticos o de diseño 


Volpi y Fiori (2012) destacan que la gráfica del periodo de retorno 
conjunto de tipo AND, mostrada más adelante como Figura 2, presenta 
una severa inconsistencia al contener, en un contexto bivariado, umbrales 
críticos univariados. Debido a lo anterior, tal gráfica se considera 
integrada por dos porciones: las dos designadas simples (naive part) y la 
correcta (proper part). Las partes rectas son las colas o rectas asíntotas 
a la parte curva. La probabilidad de ocurrencia de un evento o pareja de 
Q y V es variable en la parte curva y decrece a lo largo de la parte recta, 
aunque todos los valores definen el mismo periodo de retorno conjunto. 
En resumen, las parejas de valores de las rectas asíntotas tienen 
probabilidades de ocurrencia bajas y por ello no deben ser incluidas en 
los análisis de búsqueda de las crecientes (Q y V) críticas o severas. Para 
fines prácticos, los puntos extremos de la parte curva se pueden definir 


según su distribución empírica o cercanos al inicio de las rectas asíntotas. 


Selección de las distribuciones marginales 


El planteamiento para la selección de las distribuciones marginales fue 
muy simple y consistió en aplicar las tres FDP que se han establecido de 
referencia o de aplicación bajo precepto, que son la Log-Pearson tipo III 
(LP3), la General de Valores Extremos (GVE) y la Logística Generalizada 


(LOG). Además, se aplicaron tres distribuciones de uso generalizado: la 
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Pareto Generalizada (PAG), la Kappa y la Wakeby. Las primeras cuatro 
FDP citadas tienen tres parámetros de ajuste, y las dos últimas, cuatro y 
cinco. Con excepción de la LP3 que se aplicó con el método de momentos 
en los dominios logarítmico (WRC, 1977) y real (Bobée, 1975), el resto 


se ajustó con el método de momentos L (Hosking € Wallis, 1997). 


La selección de la FDP más conveniente se basó en el valor que cada una 
genera con una probabilidad de no excedencia del 1 %; es decir, un valor 
muy bajo que no supere a los valores mínimos observados para evitar 
probabilidades marginales (u, v) negativas. Este criterio de selección es 
el requerido en registros de Q y V que presentan valores mínimos muy 


reducidos, en comparación con sus extremos máximos. 


Además, se tomaron en cuenta a los errores estándar de ajuste (EEA) y 
absoluto medio (EAM), así como la magnitud de las predicciones en los 


periodos de retorno mayores a 500 años. 


Errores de ajuste 


El primer criterio aplicado para la selección de la mejor FDP a unos datos 
o serie disponible fueron los llamados errores de ajuste (Kite, 1977; 
Willmott € Matsuura, 2005; Chai 8 Draxler, 2014). Este criterio y el 
descrito para evitar probabilidades negativas permitirán adoptar una 
distribución adecuada entre los modelos: LP3, GVE, LOG, PAG, Kappa y 
Wakeby. 


Cambiando en las ecuaciones (28) y (29) las probabilidades observadas 


por los datos ordenados de la serie analizada (x; O y;), y las probabilidades 
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calculadas por los valores estimados con la FDP que se prueba o 
contrasta, se obtienen el error estándar de ajuste (EEA) y el error absoluto 
medio (£EAM). Los valores que se estiman (%; o 9) se buscan para la misma 
probabilidad de no excedencia, asignada a los datos por la fórmula 
empírica de Gringorten (Ecuación (26)). 


Datos por procesar 


Para realizar el contraste de las cuatro FC que tienen dependencia 
significativa en su extremo derecho (Ay), se procesaron los 16 registros 
conjuntos de gasto máximo (Q) y volumen (V) de las crecientes anuales 
que expusieron Aldama et al. (2006), los cuales suman 620 parejas de 
datos en 15 embalses importantes de México y uno en proyecto (La 


Parota). 


Además, se incluyó el registro conjunto de Q y V anuales de las crecientes 
registradas en la estación hidrométrica La Cuña, expuesto por Gómez et 
al. (2010) y Campos-Aranda (2022). También se procesó el registro de Q 
y V de las crecientes de entrada a la presa Malpaso, mostrado por 


Domínguez y Arganis (2012). 


En la Tabla 2 se expone el registro de los 26 datos conjuntos de Q y V de 
las crecientes de entrada a la presa Adolfo López Mateos en el río Humaya 
del estado de Sinaloa, México. Tal registro procede de Aldama et al. 


(2006) y se procesó como aplicación numérica. 
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Tabla 2. Gastos máximos, volúmenes y sus números de orden 
bivariados en las 26 crecientes anuales de entrada a la presa Adolfo 


López Mateos (Humaya), México (Aldama et al., 2006). 


Año e Ñ NM; Año A y NM; 
(m/s) | (Mm) (m*/s) | (Mm) 
(1974 | 600 | 140 | 4 |1987 | 401 | 160 | 6 | 
(1975 290 | 200 | 4 |1988 | 305 | 150 | 3 | 
(1976 316 | 70 | 3 |1989| 581 — 170 | 8 | 
(1977 242 | 50 | 1 |1990 | 2035 | 1150 | 22 | 
(1978 248 | 420 | 4 |1991| 488 — 190 | 9 | 
[1979 891 | 400 | 15 ¡192 | 242 | 310 | 3 | 
1980 461 190 7 | 1993 487 160 7 
1981 | 9245 4580 | 26 | 1994 933 150 6 
1982 | 2388 1410 | 23 | 1995 127 70 1 
1983 277 160 3 (1996 | 4490 1240 | 24 
1984 481 290 9 | 1997 904 160 9 
1985 614 450 16 | 1998 | 2529 1020 | 22 
1986 | 1064 550 20 | 1999 | 1193 380 7 
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Test de Wald-Wolfowitz 


Esta prueba no paramétrica ha sido utilizada por Bobée y Ashkar (1991), 
Rao y Hamed (2000), y Meylan et al. (2012) para verificar independencia 
y estacionariedad en registros de gastos máximos anuales (X;). Por lo 
anterior, se propuso aplicar esta prueba a los registros conjuntos de Q y 


V, que deben ser muestras aleatorias. 


Resultados y su discusión 


Contraste de funciones Cópula con 2; > O 
Verificación de la aleatoriedad 


En la Tabla 3 se muestran en orden progresivo de tamaño de cuenca (4 
en km2) los 16 registros conjuntos de Q y V procesados, tomados de 
Aldama et al. (2006) y, al final, los dos registros conjuntos también 
procesados. En la columna 4 se expone el número de datos (n) de cada 
registro, y en las columnas 5 y 6, el valor del estadístico U del test de 
Wald-Wolfowitz, cuyo valor absoluto debe ser menor de 1.96 para que tal 


registro sea aleatorio. 


160 


Tecnología y ciencias del agua, ISSN 2007-2422, 


Open Access bajo la licencia CC  BY-NC-SA 4.0 | 14(5), 120-188. DOI: 10.24850/j-tyca-14-05-03 


(https: //creativecommons.org/licenses/by-nc-sa/4.0/) 


OPEN Que: 


Tecnología y 


CienciaszAgua 


0) Check for updates 


Tabla 3. Datos generales y resultados del test de Wald-Wolfowitz en los 


18 registros conjuntos de Q y V procesados. 


Test WW 
No. Nombre del embalse o registro A(km?) | n 
Q v 

1 Madin 171 30 | -0.757 | -0.335 
2 Eustaquio Buelna (Guamúchil) 1630 37 | (4.156) | (2.597) 
3 Josefa Ortiz de Domínguez (El Sabino) 2 268 34 | (1.986) | 1.529 
4 Abelardo L. Rodríguez 2 430 44 | (2.220) | (2.346) 
5 Sanalona 3 250 38 | 0.337 | -1.555 
6 Ignacio Allende (La Begoña) 4 984 27 | 0.287 1.241 
7 La Parota 7 067 31 | -0.760 | -0.532 
8 Benito Juárez (El Marqués) 9 697 26 | -0.601 | 0.427 
9 Adolfo Ruiz Cortines (Mocúzari) 10 719 61 0.147 1.189 
10 | Adolfo López Mateos (Humaya) 10 972 26 0.677 1.252 
11 | Belisario Domínguez (La Angostura) 18 537 33 | 0.955 0.247 
12 | La Boquilla 21 003 64 | -0.084 | -0.251 
13 | Luis Donaldo Colosio (Huites) 26 020 52 | -0.091 | -0.602 
14 | Venustiano Carranza (Don Martín) 31 034 52 0.726 1.054 
15 | Plutarco Elías Calles (E/ Novillo) 58 280 22 0.637 0.496 
16 | Adolfo López Mateos (Infiernillo) 108 000 | 43 | -0.707 | -0.022 
17 | Estación hidrométrica La Cuña 19097 | 55 | 0.284 0.213 
18 | Presa Malpaso (Netzahualcóyotl) 34 800 47 | -0.666 | 0.558 
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Se deduce de los valores citados que dos registros conjuntos no son 
aleatorios: los de las presas Guamúchil y Abelardo L. Rodríguez. A tales 
registros se les aplicaron pruebas estadísticas específicas y se encontró 
en la presa Guamúchil que las series de Q y V muestran persistencia y 
exceso de variabilidad por sus valores extremos. Además, la serie de Q 
muestra tendencia ascendente. En la presa Abelardo L. Rodriguez, ambas 
series de Q y V muestran persistencia, exceso de variabilidad y tendencia 
negativa. El registro de Q de la presa El Sabino puede ser aceptado 
aleatorio, ya que sólo presenta persistencia ligera. Por lo anterior, los 
registros conjuntos de las presas Guamúchil y Abelardo L. Rodríguez se 


eliminaron de los análisis posteriores. 


Valores de las medidas de asociación 


En la columna 2 de la Tabla 4 se muestran los valores del cociente tau de 
Kendall obtenidos con la Ecuación (21). Respecto al coeficiente rho de 
Spearman, mostrado en la columna 2 de la Tabla 5, los valores que no 
tienen paréntesis resultaron iguales al calcularse con las ecuaciones (23) 
y (24). En cambio, las magnitudes entre paréntesis provienen de la 


Ecuación (23), que conduce a valores ligeramente mayores. 
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Tabla 4. Valores del cociente tau de Kendall de la dependencia 


observada (Ay “), y de las Cópulas aplicadas de Joe y de Clayton 


Asociada al registro conjunto indicado. 


Embalse o en Cópula Joe Clayton Asoc. 
Ta A 
registro E 0 ww 0 du 
Madin 0.6598 | 0.7356 | 4.687 | 0.8406 | 3.8784 0.8363 


El Sabino 0.5936 | 0.6776 | 3.748 | 0.7969 | 2.9211 0.7888 


Sanalona 0.5733 | 0.6053 | 3.525 | 0.7827 | 2.6867 0.7726 


La Begoña 0.6866 | 0.7831 | 5.177 | 0.8567 | 4.3818 | (0.8537) 


La Parota 0.5226 | 0.6724 | 3.039 | 0.7438 | 2.1892 | (0.7286) 


El Marqués [0.5385 | 0.6462 | 3.179 | 0.7564 | 2.3333 0.7430 


Mocúzari 0.7093 | 0.8024 | 5.630 | 0.8690 | 4.8797 | (0.8676) 


Humaya 0.4523 | 0.6834 | 2.522 | 0.6837 | 1.6517 0.6573 


La Angostura |0.4318 | 0.4730 | 2.397 | 0.6647 | 1.5200 0.6338 


La Boquilla 0.6091 | 0.6440 | 3.938 | 0.8075 | 3.1168 | 0.8006 


Huites 0.4495 | 0.5655 | 2.504 | 0.6811 | 1.6329 0.6541 


Don Martín 0.6456 | 0.7446 | 4.455 | 0.8317 | 3.6426 0.8267 


El Novillo 0.4286 | 0.6166 | 2.378 | 0.6616 | 1.5000 0.6300 


Infiernillo 0.4064 | 0.4714 ' 2.253 | 0.6398 | 1.3694 0.6028 


La Cuña 0.7199 | 0.7819 | 5.887 | 0.8750 | 5.1394 | 0.8738 


Malpaso 0.3747 | 0.4159 | 2.092 | 0.6072 | 1.1982 0.5608 
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Tabla 5. Valores del coeficiente rho de Spearman, de la dependencia 


observada (Aj ”), y de las Cópulas de Gumbel-Hougaard y t de Student 


en el registro conjunto indicado. 


Embalse o Gumbel-Hougaard t de Student 
Pn Ae 
registro 0 du du(v=2) | h(v = 10) 
Madin (0.8603) | 0.7356 | 2.9392 0.7340 0.6678 | 0.3834 
El Sabino 0.6883 0.6776 | 2.4605 0.6746 0.6072 0.2965 
Sanalona 0.7595 0.6053 | 2.3433 0.6558 0.5889 0.2725 


La Begoña 0.8107 0.7831 | 3.1909 | (0.7574) | 0.6929 | 0.4225 


La Parota 0.6968 0.6724 | 2.0946 | (0.6077) | 0.5443 | 0.2184 


El Marqués 0.7648 0.6462 | 2.1667 0.6230 0.5581 | 0.2345 


Mocúzari (0.8827) | 0.8024 | 3.4399 | (0.7768) | 0.7142 | 0.4570 


Humaya (0.7258) | 0.6834 | 1.8258 | 0.5383 | 0.4849 | 0.1563 


La Angostura 0.5836 0.4730 | 1.7600 | 0.5173 | 0.4681 | 0.1409 


La Boquilla 0.7922 0.6440 | 2.5584 | 0.6888 | 0.6212 | 0.3157 


Huites 0.5980 0.5655 | 1.8164 | 0.5354 | 0.4825 | 0.1541 


Don Martín (0.8303) | 0.7446 | 2.8213 0.7215 0.6546 | 0.3636 


El Novillo 0.5788 0.6166 | 1.7500 0.5140 0.4654 | 0.1385 


Infiernillo (0.5707) | 0.4714 | 1.6847 0.4910 0.4476 | 0.1233 


La Cuña (0.9088) | 0.7819 | 3.5697 0.7857 0.7243 | 0.4737 


Malpaso (0.5510) | 0.4159 | 1.5991 0.4574 | 0.4226 | 0.1038 
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Contraste de dependencias (1,) 


En las columnas 3 de la Tabla 4 y Tabla 5 se exponen los valores de la 
dependencia observada (A5"“) estimada con la Ecuación (36). El criterio 
de selección de la FC fue que su dependencia significativa (1,4) debería ser 
algo mayor o escasamente menor a la dependencia observada (A£"“). Con 
tal criterio, en la Tabla 4 se definieron las FC de Joe para el registro de la 
presa Adolfo López Mateos (Humaya) y la FC de Clayton Asociada para el 


registro de la presa Plutarco Elías Calles (El Novillo). 


Por otra parte, en la Tabla 5, con el criterio de selección descrito, se 
definen 10 registros, cuya dependencia en su extremo derecho puede ser 
presentada por la FC de Gumbel-Hougaard (GH), y tres que están 
comprendidos entre las FC de Clayton Asociada y de GH, cuyos valores 
de la dependencia (A,,) se han mostrado entre paréntesis. La FC de la t de 
Student reproduce la dependencia del registro de la presa Malpaso, con 


v=2. 


Búsqueda de las distribuciones marginales 
Distribución de los gastos máximos anuales 


En la Tabla 6 se exponen los errores de ajuste y predicciones (m?/s) 
obtenidas con las seis distribuciones aplicadas al registro de gastos 
máximos de la Tabla 2, cuyo valor mínimo registrado es 127 m3/s. De 


estas seis FDP, las únicas que definen valores con probabilidad de no 
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excedencia del 1 % menores al gasto citado son la GVE y la LOG. Tomando 


en cuenta los errores de ajuste, se adopta la GVE. 


Tabla 6. Errores de ajuste y predicciones (m3/s) de las seis 


distribuciones aplicadas en el registro de gastos máximos anuales de las 


crecientes de entrada a la presa Adolfo López Mateos (Humaya), México. 


Periodos de retorno en años 


ab EEA EAM 

(m3/s)  (m?/s) 50 100 500 1 000 5000 | 10000 
LP3 230.8 132.6 | 9062 | 15192 | 48746 | 79778 | 247 228 | 400 939 
GVE 675.5 279.8 6663 | 10381 | 28653 | 44 244 | 121000 | 186 452 
LOG 745.8 301.9 6364 | 9961 28041 | 43766 | 123060 | 192 052 
PAG 631.9 255.6 6824 | 10412 | 26828 | 40002 | 100 278 | 148 626 
KAP 585.1 243.0 7072 | 10697 | 26739 | 39253 | 94582 | 137 656 
WAK | 676.7 286.5 6770 | 10379 | 27111 | 40694 | 103 672 | 154 793 


Los parámetros de ubicación (u:), escala (ai) y forma (k1) de la 
distribución GVE adoptada son 447.7865, 372.8246 y -0.6236228, 


respectivamente, cuya ecuación es: 


Open 
(https: //creativecommons.org/licenses/by-nc-sa/4.0/) 


Access 
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licencia 26€ 


BY-NC-SA 
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Distribución de los volúmenes anuales 


En la Tabla 7, similar a la anterior para los volúmenes anuales de la Tabla 


2, cuyo valor mínimo es 50 Mm3, nuevamente las únicas FDP que no 


generan probabilidades marginales negativas son la GVE y LOG. Se 


adopta la GVE por mostrar menores errores de ajuste. 


Tabla 7. Errores de ajuste y predicciones (m3/s) de las seis 


distribuciones aplicadas al registro de volúmenes anuales de las 


crecientes de entrada a la presa Adolfo López Mateos (Humaya), México. 


EEA EAM Periodos de retorno en años 

FDP 

(m?/s) | (m*/s) [50 100 500 1000 | 5000 | 10000 
LP3 189.8 96.6 4061 | 6726 | 20745 | 33252 | 97487 | 153 958 
GVE 370.6 138.8 | 3021 | 4726 | 13163 | 20401 | 56 269 | 87 019 
LOG 400.4 144.6 2881 | 4528 | 12856 | 20134 | 57053 | 89 331 
PAG 355.1 133.1 3091 | 4739 | 12345 | 18494 | 46869 | 69 796 
KAP 338.6 135.1 3207 | 4875 | 12327 | 18186 | 44336 | 64 877 
WAK 375.3 146.4 3075 | 4 729 12 428 18697 | 47 873 71 630 


Los parámetros de ubicación 


(u2), escala (a2) y forma (k>2) de la 


distribución GVE adoptada son 193.4795, 167.2692 y -0.6288097, 


respectivamente, cuya expresión es: 


Open 


Access bajo 


la licencia 
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1, 1/42 
F(y) = exp (- [1 — k, al ] (68) 


Selección y ratificación de la FC 


El procesamiento bivariado de los datos de la Tabla 2 condujo a los 
siguientes tres indicadores de asociación: rxy = 0.9670, rt, = 0.4523 y Pn 
= 0.6738. La Ecuación (22), con n = 26 y el tau citado origina un valor 
de 3.24; por lo cual, el cociente de Kendall es significativo. La Ecuación 
(25) genera un valor de 2.26 y entonces el coeficiente de Spearman 


también es significativo. 


Por otra parte, en la Tabla 8 se muestran los indicadores estadísticos de 
ajuste que se obtuvieron al aplicar las FC de Clayton, Frank, Plackett, t- 
Student, Gumbel-Hougaard (GH), Clayton Asociada y Joe. En las 
ecuaciones (28) a (30), las probabilidades bivariadas empíricas se 
estimaron con la Ecuación (27) y las teóricas con las ecuaciones (2), (4), 
(9), (41), (7), (18) y (37). Los valores reportados para la Cópula t- 


Student corresponden a un v = 6, donde ya estabilizaron los errores de 


TI Tn 


ajuste con O = sen ( ) = 0.6522, según la Ecuación (42). 


2 
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Tabla 8. Indicadores estadísticos del ajuste de las funciones Cópula 
indicadas en las crecientes anuales de entrada a la presa Adolfo López 


Mateos (Humaya), México. 


, Núm. | Núm. 
Cópula 0 EME EAM MDP MDN Ay 
DP DN 
Clayton (C.) 1.6517 | 0.0529 | 0.0365 14 1.2 0.1505 -0.0779 0.0000 
Frank 4,9310 | 0.0461 | 0.0321 15 di 0.1334 -0.0701 0.0000 
Plackett 11.079 | 0.0441 | 0.0315 11 15 0.1291 -0.0727 0.0000 
t-Student 0.6522 | 0.0482 | 0.0347 13 13 0.1415 -0.0643 0.2641 
GH 1.8258 | 0.0473 | 0.0354 14 1.2 0.1352 -0.0601 0.5383 
C. Asociada 1.6517 | 0.0475 | 0.0362 14 1.2 0.1277 -0.0549 0.6573 
Joe 215220 10.0477 7 10105865 15 pel 0.1267 -0.0533 | 0.6837 
Significado de los nuevos acrónimos: 
DP, DN = diferencias positivas y negativas. 
MDP, MDN = máxima diferencia positiva y negativa. 
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Con base en los resultados de la Tabla 8, no existe dificultad para 
seleccionar como mejor FC la de Plackett, para los datos de la Tabla 2, 
debido a que reporta los indicadores estadísticos (EME y EAM) más 


reducidos (mostrados en cursivas). 


Sin embargo, al tomar en cuenta que la dependencia observada (A$5"“) 
alcanzó un valor de 0.6834, la FC que se debe seleccionar es la Joe, como 
ya se expuso en la Tabla 4. En la Tabla 8 se observa que el ajuste de la 
FC de Joe a los datos es mejor que el peor de todos, que lo definió FC de 
Clayton. Es destacable que la FC de Joe origina los valores más bajos de 
las diferencias máximas entre probabilidades empíricas y teóricas 


(mostradas en cursivas y sombreadas). 


En la Tabla 9 se muestran las probabilidades de no excedencia bivariadas, 
empíricas observadas (w?) y teóricas calculadas (wf) con la FC de Joe. 
También se indican sombreadas las diferencias máximas positiva y 


negativa. 
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Tabla 9. Probabilidades de no excedencia conjuntas y sus diferencias 


calculadas con la FC de Joe para las crecientes de anuales de entrada a 


la presa Adolfo López Mateos (Humaya), México. 


Núm wi wi Diferencias | Núm. wi wi Diferencias 
1 0.1363 | 0.1896 -0.0533 14 |0.2129 | 0.1619 0.0509 
2 0.1363 | 0.1300 0.0063 15 |0.0980 | 0.1051 -0.0071 
3 0.0980 | 0.0311 0.0669 16 |0.2894 | 0.2392 0.0502 
4 0.0214 | 0.0101 0.0114 17 0.8254 | 0.8643 -0.0389 
5 0.1363 | 0.1378 -0.0015 18 |0.3277 | 0.2369 0.0909 
6 0.5574 | 0.5663 -0.0089 19 |0.0980 | 0.1215 -0.0234 
7 0.2511 | 0.2249 0.0263 20 0.2511 | 0.1955 0.0556 
8 0.9786 | 0.9851 -0.0066 21 |0.2129| 0.2460 -0.0331 
9 0.8637 | 0.8939 -0.0302 22 |0.0214| 0.0051 0.0163 
10 | 0.0980 | 0.0973 0.0007 23 |0.9020| 0.9195 -0.0175 
11 0.3277 | 0.3210 0.0067 24 0.3277 | 0.2666 0.0612 
12 0.5957 | 0.4690 0.1267 25 0.8254 | 0.8767 -0.0513 
13 | 0.7489 | 0.6690 0.0798 26 0.6340 | 0.6058 0.0282 


Por otra parte, la Ecuación (52) define a D, = 0.2663 y como la diferencia 


máxima absoluta de la Tabla 9 es 0.1267, el test de Kolmogorov-Smirnov 


ratifica a la FC de Joe adoptada. El coeficiente de correlación (rxy) entre 


las probabilidades empíricas y las teóricas estimadas con la FC de Joe 
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resultó de 0.9897; por lo tanto, aceptable. El contraste gráfico entre 
ambas probabilidades para ratificar su adopción se expone en la Figura 1 


para los datos de la Tabla 9. 


órica 


r 


Probabilidad de no excedencia te 


0.0 0.2 0.4 0.6 0.8 1.0 
Probabilidad de no excedencia empírica 


Figura 1. Contraste gráfico de probabilidades conjuntas (Q, V) 
calculadas con la FC de Joe para las crecientes anuales de entrada a la 


presa Adolfo López Mateos (Humaya), México. 
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Gráficas del periodo de retorno bivariado T,, 


Los periodos de retorno bivariados de tipo AND se estiman con la Ecuación 
(54). Para sus estimaciones con valores del Tyy de 500, 1 000, 5 000 y 10 
000 años se aplicó la FC de Joe. De manera arbitraria se seleccionan 
volúmenes y gastos pico para obtener sus probabilidades de no 
excedencia marginales (ecuaciones (67) y (68)) y conjunta (Ecuación 
(54)). En la Tabla 10 se muestran resultados para definir las cuatro 
gráficas de la Figura 2. 


Tabla 10. Parejas de gasto pico y volumen anual utilizadas para definir 
las gráficas del periodo de retorno conjunto tipo AND con la FC de Joe 
en las crecientes de entrada a la presa Adolfo López Mateos (Humaya), 
México. 
Txy Txy Txy Txy 


500 años 1 000 años 5 000 años 10 000 años 


VMm3 |Qm?/s | VMm? | Q m?/s | VMm* Q m?/s VMm3 | Qm?/s 


0 28 653 0 44 244 0 121 000 0 186 452 


1 000 28 636 2 000 44 200 4 000 120 940 5 000 186 420 


2 000 28 571 4 000 44 020 8 000 120 735 10000 | 186 300 


3 000 28 439 6 000 43642 | 12 000 120 290 15000 | 185 840 


5 000 27 903 8 000 43002 | 16 000 119 500 20000 | 185 250 


6 000 27453 | 10000 | 42014 | 20 000 118 390 25000 | 184 150 


7 000 26836 | 12000 | 40542 | 24 000 116 830 30000 | 182 800 
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Txy Txy Txy Txy 
500 años 1 000 años 5 000 años 10 000 años 


VMm3 |Qm?/s | VMm? | Q m?3/s | VMm*3 Q m?/s VMm3 | Qm?/s 


8 000 26000 | 13000 | 39560 | 28 000 114 700 35000 | 181 120 


9 000 24 858 | 14000 | 38 348 | 32 000 111 900 40000 | 178 750 


10000 | 23266 | 15000 | 36850 | 36 000 108 140 45000 | 175 800 


10500 | 22222 | 16000 | 34960 | 40 000 103 060 50000 | 172 050 


11000 | 20940 | 17000 | 32520 | 42 000 99 850 52 500 | 169 960 


11500 | 19315 | 18000 | 29212 | 44 000 96 060 55000 | 167 520 


12000 | 17165 | 18500 | 27050 | 46 000 91 430 57 500 | 164 740 


12500 | 14046 | 19000 | 24 335 | 47 000 88 770 60000 | 161 570 


12750 | 11745 | 19500 | 20720 | 48 000 85 800 62500 | 157 970 


13 000 8 125 19 750 | 18330 | 49 000 82 460 65000 | 153850 


13 050 6 990 20000 | 15170 | 50 000 78 670 67500 | 149 310 


13 100 5 490 20100 | 13530 | 51 000 74 250 70000 | 143810 


13 163 0 20200 | 11520 | 52 000 69 140 72000 | 138820 


20 300 8 730 53 000 62 840 74000 | 133200 


20 350 6 450 53 500 59 250 76000 | 126 530 


20 401 0 54 000 54 985 78000 | 118 720 


54 500 50 160 80000 | 109 000 


55 000 44 100 81000 | 103 580 


55 500 35 840 82 000 96 950 
56 000 23 600 83 000 89 650 


56 269 0 84 000 80 650 
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Txy Txy Txy Txy 
500 años 1 000 años 5 000 años 10 000 años 
VMm3 |Qm?/s | VMm? | Q m?3/s | VMm* Q m?/s VMm3 | Qm?/s 


85 000 68 990 


86 000 53 500 


87 019 0 


100 


90 A7'(Q ,V )=10,000 años 


a 


Volumen anual (10? Mm”) 


60 AT" (Q ,V )=5,000 años 


20 40 60 80 100 120 140 160 180 200 


Gasto máximo anual (10? m*/s) 


Figura 2. Gráficas de los cuatro periodos de retorno conjunto Ty, de 


diseño obtenidas con la FC de Joe en las crecientes anuales de entrada a 


la Presa Adolfo López Mateos (Humaya), México. 
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Selección de los eventos de diseño 


En la Figura 2 o en la Tabla 10 se pueden seleccionar infinitas parejas de 
Q y V, que satisfacen el periodo de retorno conjunto de diseño y que se 
definen como subgrupo de parejas críticas por estar dentro de la porción 
curva de cada gráfica del Ty, fuera de las rectas asíntotas (Volpi 8 Fiori, 
2012). 


Las combinaciones de gasto pico y volumen que tienen el mismo periodo 
de retorno bivariado de diseño establecen crecientes o hidrogramas que 
producirán diferentes efectos en el embalse que se diseña o revisa, 
adoptando por seguridad, el que genera las condiciones más críticas. Para 
formar cada hidrograma de diseño existen métodos teóricos y empíricos 
(Aldama, 2000; Aldama et a/., 2006; Campos-Aranda, 2008; Gómez et 
al., 2010; Gráler et a/., 2013). 


Contraste de los periodos de retorno conjuntos 


En la Tabla 11 y Tabla 12 se muestran los cálculos realizados para llevar 
a cabo la verificación de las ecuaciones (55) y (62) con base en las FC de 
Gumbel-Hougaard y de Joe. Se observa y comprueba en todos los casos 
que ambas FC cumplen con la Ecuación (55). Para el periodo de retorno 
secundario, la Ecuación (62) se cumple exclusivamente con la FC de Joe 


en los tres últimos periodos de retorno analizados. 
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Tabla 11. Periodos de retorno univariados y bivariados estimados con la 
FC de Gumbel-Hougaard en las crecientes de entrada a la presa Adolfo 


López Mateos (Humaya), México. 


FC de Gumbel-Hougaard Txy secundario 
Tx y Ty 
CIFXOO,Fr(y)] Txy Txy Kc(s) g 
50 0.9709003 34.4 91.7 0.990844 109.2 
100 0.9854162 68.6 184.6 0.995450 219.8 
500 0.9970778 342.2 927.9 0.999094 | 1 104.2 
1000 0.9985386 684.3 | 1856.7 | 0.999547 | 2 209.6 
5 000 0.9997076 3 420.4| 9294.9 | 0.999910 |11 059.5 
10 000 0.9998538 6 836.7 |18 600.0 | 0.999955 | 22 104.4 


Tabla 12. Periodos de retorno univariados y bivariados estimados con la 
FC de Joe en las crecientes de entrada a la presa Adolfo López Mateos 


(Humaya), México. 


FC de Joe T secundario 
Tx y Ty 
CIFXCO,Fr(y)] Txy Txy Kc(s) g 
50 0.9736739 38.0 73.1 0.987932 82.9 
100 0.9868368 76.0 146.3 0.993994 | 166.5 
500 0.9973674 379.9 731.4 0.999022 |1 022.7 
1000 0.9986837 759.7 | 1462.7 | 0.999000 |1 000.0 
5000 | 0.9997368 |3799.2| 7313.5 | 0.999800 |5000.7 | 
10000 | 0.9998683  |7594.9|14627.0| 0.999900 |9 998.3 | 
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Conclusiones 


En México, los seis trabajos citados: Aldama (2000), Ramírez-Orozco y 
Aldama (2000), Escalante-Sandoval y Reyes-Chávez (2002), Aldama et 
al. (2006), Gómez et al. (2010), y Domínguez y Arganis (2012) ponen de 
manifiesto los esfuerzos realizados para desarrollar métodos confiables 
de estimación de crecientes, que brinden seguridad a las obras hidráulicas 
rurales y urbanas que son dimensionadas y/o revisadas con ellas, al 
buscar estimar crecientes bivariadas lo más aproximadas posible; este 


estudio se suma a los esfuerzos citados. 


La ventaja fundamental de usar las funciones Cópula (FC) en el análisis 
de frecuencias de crecientes bivariado (AFCb) radica en construir 
fácilmente la distribución de probabilidades conjunta con base en las 
distribuciones univariadas marginales iguales o diferentes adoptadas, 
previa estimación de la dependencia entre las variables aleatorias: gasto 
máximo (Q) y volumen escurrido (VW) de las crecientes anuales 


disponibles. 


La aplicación de las FC en los AFCb ha mostrado que la selección de la 
Cópula idónea se basa en reproducir la dependencia del extremo derecho 
observada, lo cual garantiza una estimación más aproximada de las 
predicciones asociadas con altos periodos de retorno conjuntos. Por lo 
anterior, ahora se prueban FC que exhiben, de más a menos, varios 
grados de dependencia en su cola derecha, por las familias de las Cópulas 


de Joe, Clayton Asociada, Gumbel-Hougaard y t de Student. 
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El contraste realizado con 18 registros anuales de Q y V recopilados 
detectó dos no aleatorios, que fueron eliminados. Para el resto, la FC de 
Joe se sugiere en un registro procesado como aplicación numérica. La FC 
de Clayton Asociada se adopta en otro registro y la FC de Gumbel- 
Hougaard se propone en 10 registros. En tres registros su dependencia 
observada está comprendida entre la que exhiben estas dos FC citadas. 
Por último, se encontró que la FC de la t de Student reproduce la 


dependencia del registro de la Presa Malpaso, con v = 2 (ver Tabla 5). 


Las FC que no muestran tal dependencia, como las Clayton, Frank y 
Plackett, también se aplican a los datos disponibles de Q y V, y permiten 
validar la calidad del ajuste logrado por la FC probada. Lo anterior está 
respaldado por el hecho de que todas las FC citadas y expuestas muestran 


un ajuste similar a los datos (ver Tabla 8). 


El procesamiento de los 26 datos anuales de Q y V de las crecientes de 
entrada a la presa Adolfo López Mateos (Humaya) del estado de Sinaloa, 
México, permitió mostrar el desarrollo numérico de la teoría expuesta y 
describir con detalle las Figuras 1 y 2, que engloban la aceptación de la 
FC adoptada y la gráfica básica de los resultados o predicciones de Q y V 


del periodo de retorno conjunto tipo AND. 
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de crecientes. Sus comentarios permitieron destacar los esfuerzos 
realizados en México para contar con métodos confiables de estimación 


de crecientes de diseño en grandes embalses. 
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